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Abstract 

Several weeks after the explosion of supernova (SN) 1987A, the UV flash 
of the SN illuminated a ring-like structure in the circumstellar material at 
about 0.65 ly from the SN. The interaction between the stellar winds from 
the SN progenitor is considered to be the candidate for the formation of the 
circumstellar structure. In the case that the stellar winds are spherically 
symmetric, the interaction should result in a shell- like structure. How- 
ever, in this paper we show that the magnetic field in the winds causes an 
anisotropy which leads to the formation of a ring-like structure. When 
the fast wind of the blue supergiant phase of the progenitor sweeps up the 
surrounding slow wind of the red-supergiant phase, the magnetic field as 
well as the wind material are piled up in the interaction region. Since the 
magnetic energy increases in proportion to the square of the amplitude, the 
magnetic field exhibits its effect prominently at the interaction region; due 
to the magnetic pressure force the material at lower latitudes is compressed 
into a ring-like structure. It is suggested that this magnetic process can 
also explain the newly observed pair of rings of the SN 1987A nebula. 
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1. Introduction 

It is widely believed that the progenitor of SN 1987A was a red super- 
giant (RSG) which then evolved to a blue supergiant (BSG) just prior to 
the explosion. In such a case, the circumstellar material surrounding the 
SN would comprise two types of stellar winds: a dense and slow RSG gen- 
erated wind far removed from the SN, and a dilute, fast BSG wind near to 
it. When both winds are isotropic, as the BSG wind expands hydrodynam- 
ically into the RSG wind, a shell-like density enhancement is expected to 
form in the interaction region between the two winds. Actual observations, 
however, do not support this simple picture. Results from the European 
Southern Observatory (ESO, Wamplcr et al. 1990) and the Hubble Space 
Telescope (HST, Jakobson et al. 1991 and Panagia et al. 1991) reveal an 
unexpected ring-like structure, rather than a shell-like structure, at about 
0.65 light-year (ly) from SN 1987A. 

In the early literature, ring-like structures were discussed in terms of 
planetary nebulae. Moreover, such nebulae are also believed to generate 
a fast wind propagating into a slower wind. Several workers [e.g., Kwok 
(1982), Kahn and West (1985), Balick, Preston, and Icke(1987), Socker and 
Livio(1989)] have used hydrodynamic (HD) schemes to study the interaction 
region along with the assumption of a relatively high-density concentration 
along the equator of the slow wind. They found that the interaction between 
the slow and fast winds enhances the density-anisotropy in the slow wind in 
such a way that a ring-like structure is formed. Based on this idea, Luo and 
McCray (1991) suggest that the ring-like structure of the SN 1987A nebula 
resulted from an anisotropy ratio of equatorial-to-polar density of 5 in the 
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RSG wind (slow wind) before the interaction. Wang and Mazzali (1992) 
found that the outward propagating shell from the fast-slow boundary can 
encounter and substantially compress the anisotropy of a weak (ratio 1.2) 
slow-wind. 

All of the above HD studies, however, seem to have neglected the pressure 
of the compressed gas at the interaction region, which would prevent the 
concentration of gas. If this effect is taken into account, a much higher 
anisotropy ratio would be necessary for ring formation. In fact, a recent 
full-HD computer simulation by Blondin and Lundqvist (1993) shows that 
a very high-density anisotropy of the order of 20 in the background RSG 
wind is required for ring formation. In addition, a new image taken by 
the Hubble Space Telescope shows an additional pair of rings, which has a 
co-axis with the original ring (McCray, Lin 1994). This new image seems 
to be inaccessible in terms of the HD point of view. 

In this paper we propose an MHD model in which the magnetic pressure 
effect can account for the formation of the ring-like structure. The main 

purpose is to explain the original ring, but as will be shown below, a pair 
of rings is also reproduced in our MHD simulation, though the geometrical 
scale should be adjusted by a further fine tuning of the model. 

Based on an evolutionary scenario of the progenitor's magnetic field and 
rotation, the RSG wind is likely to posess such a magnetic field strength 
and rotation rate that the magnetic pressure causes a density anisotropy. 
It is shown for a magnetic stellar wind from a dipole source field (Washimi, 
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Shibata 1993) that the circumstellar poloidal magnetic field {Bp, the field on 
the meridional plane) configuration becomes a split monopole configuration; 
the field is oriented to the radial direction, but its polarity changes in the 
equatorial region, and, thus, a magnetic neutral sheet is formed along the 
equator. The toroidal magnetic field {B^, the field of the toroidal compo- 
nent around the stellar rotation axis, or in other words, the field vertical 
to the meridional plane) is proportional to ^ ■ Bp ■ R ■ sin 0, where O is 
the stellar rotation speed, R the distance from the star, and the colati- 
tude angle. Since Bp and B^p decrease as R~^ and R^^, respectively, B^p is 
much stronger than Bp in circumstellar space far from the stellar surface. 
Thus, a thin magnetic neutral sheet of the toroidal field is formed around 
the equator where the density is enhanced to maintain the total pressure 
balance, as shown by Washimi and Shibata. During the dynamical stage 
of the sweeping up of the RSG wind by the BSG wind, the toroidal field is 
piled up in the interaction region. Since the energy increases as the square 
of the intensity, this piling up leads to a drastic increase in the magnetic 
energy. Due to all of the above reasons, the magnetic pressure force of 
the toroidal field compresses the gas (plasma) at the equatorial interaction 
region and maintains a localized high-density distribution over a long-time 
period by pressure balancing the ring of the enhanced gas. 

2. Basic Model and Simulation Method 

We perform an MHD computer analysis in order to confirm this physical 
process. Both the RSG and BSG winds are assumed to be obeyed by the 
ideal MHD system of equations, i.e.. 
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dp/dt + V • (pv) = 0, 

d^f/dt + (v ■ V)v = -(l/p)VP + (l/p)J X B 
dP/dt + (v V)P = -7PV • V, 
d'B/dt = V X (v X B), 



(1) 
(2) 
(3) 
(4) 



and 



/xqJ = V X B, 



(5) 



where ijlq is the permeabihty and the polytropic index 7 is 5/3. Here, the 
cyhndrical coordinates (r, cf), z) with the z-axis ahgned along the rotation 
axis of the progenitor, are used and axial symmetry is assumed. The flow 
velocity, v, the magnetic field, B, and the current, J, in equations (l)-(5) 
are three-dimensional vectors. 

Once a magnetic field is introduced, we may have a wide variety of parame- 
ters for the wind structure, the magnetic field configurations, and the latitu- 
dinal dependence of the wind velocity and density. In the present analysis, 
we attempt to show the ring formation using the very simple, but reasonable, 
parameters, as follows. The mass-loss rate, speed, and the temperature are 
10"^ MQ^/r-^ V^^'^ = 10 km s'^ and 5.0 x 10^ K (at 0.1 ly), respectively 
for the RSG wind, and are 3.0 x lO-^M© yr''^, V'^^'^ = 550 km s~'^ and 
1.3 X 10^ K (at 0.1 ly), respectively, for the BSG wind. Except for the 
temperature of the BSG wind, the other parameters are similar to those of 
Shigeyama and Nomoto (1990) and Luo and McCray (1991). The relatively 
high temperature of the BSG wind derives from a technical convenience to 
sustain the Rankinc-Hugoniot shock condition at the reverse shock in our 
computation. Since the ram pressure is much higher than the thermal pres- 
sure in the case of BSG wind, this choice does not lead to any unphysical 
results. The radius and the rotation speed of the progenitor are chosen to 



be R§^^ = 3.1 X 10^1 m and n^^^ = 3.5 x lO'^ rad s'^, for the RSG star; 
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whereas R^'^^ = 4.9 x 10^° m and n^^^ = 3.5 x 10"^ rad s'^ for the BSG 
star. 

Indeed, a reahstic estimate of the angular velocity as well as the magnetic 
field strength of the progenitor is very important for ring formation. There- 
fore, we need a detailed consideration of the stellar evolution for choosing 
these parameters. In general, main-sequence stars bluer than the spec- 
tral type F5 are known as rapid rotators (typically Veq ~ 100 km s'^), 
while redder stars are slow rotators (Tassoul 1978). This change in rota- 
tion along with spectral types is also observed for sub-giants and giants, 
where the change is seen around GO type (for supergiants, it is obscured 
by the convection motion). This fact suggests that the on-set of surface 
convection triggers the dynamo, and then a strong magnetic breaking of the 
stellar rotation. We can thus infer a dynamo-onset-line on the HR diagram. 

As a main-sequence star, the progenitor would rotate at typically Rq^ ~ 
100 km s-^ {n ~ 2.5 x 10"^ rad and i?o ~ 4 x 10*^ m), where Ro 
is the radius of the progenitor, which changes with stellar evolution. As 
the progenitor evolved towards RSG, it span down due to an increase in 
the moment of inertia. Taking the moment of inertia given by Eriguchi 
et al. (1992) and assuming the conservation of angular momentum, we 
have RqQ ~ 3 km, (J7 ~ 10~® rad s~^), which is consistent with the 
observed rotation velocity of typical supergiants (Gray 1991). Accord- 
ing to the suggested evolution (Saio et al. 1988), the progenitor passed 
the dynamo-onset-line just at its RSG phase, and would suffer magnetic 
breaking. If AM is the lost mass after the development of the dynamo, 
the corresponding loss of angular momentum is given by AL = AMCIR\, 
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where Ra is the averaged radius of the Alfven surface. The angular mo- 
mentum of the star is denoted by L = JMRq^}, where / ~ 4.6 x 10~^ 
(Eriguchi et al. 1992). Therefore, a prominent spin-down takes place when 
AL/L = {AM/0.1MQ){RA/2.7Rof ~ 1, i.e., when the magnetic field is 
amplified so that Ra increases up to ~ 2.7 Rq. At this stage, the corotating 
velocity Ra^ becomes ~ 8 km s~^, which is comparable with the expected 
velocity of the RSG wind. Therefore the centrifugal driving force is as 
important as other driving forces, such as due to the Alfven wave. An im- 
mediate corollary of this fact is that the magnetic pressure B^/2jj,Q beyond 
the Alfven surface is comparable with the dynamic pressure. 

Thus, it is very likely that the magnetic pressure produces a density 
anisotropy in the RSG wind. For a numerical simulation, we regard a sug- 
gestion based on the synchrotron radio emmission of SNR 1987A, and that 
the ring material was emitted at the later phase of RSG; we thus chose a set 
Bp = 12{Ro/R)'^ Gauss and O = 3.5 x 10~^ rad s~^, for which the ratio of 
the magnetic energy to the flow energy is ^ 0.9. By using this set of Bp and 
Q, we obtain a field strength on the order of 10~^ Gauss in the shocked gas, 
which is suggested by the synchrotron emission (Staveley-Smith ct al. 1992). 
Another crude estimate of the field strength on the RSG star surface may 
be provided by equipartition. According to the projenitor model by Saio 
et al. (1988), three of the radiation-energy densities, the convectional one 
and the thermal one of the gas, are comparable to each other. If the mag- 
netic energy is equated with these densities {B^/Stt ^ econv ~ o-T^ ~ nksT) 
we have 10 Gauss. However, any other set of the field strength and 
the angular velocity with the similar energy ratio could be taken for the 
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ring formation. The important point is that the condition for the density 
anisotropy is equivalent to the condition for the magneto-centrifugal driven 
wind, which is very likely to have occurred in the RSG phase. 

The poloidal magnetic field of both the RSG and BSC winds is assumed 
to be a split monopole followed by the previous discussion; the field is 
radially oriented and the colatitudinal dependence of the field is assumed to 
be Bp = bp/B? for < 9 < Oq, where hp is a constant, with a decrease to zero 
on the equator, so that Bp = bp ■ {9 ~ 90°)/(6lo - 90°)//?^ for 6*0 < ^ < 90°, 
and changes its sign in the southern-hemisphere. Wc assume = 85° 
in this simulation. The BSG poloidal field is prescribed by the magnetic 
flux conservation: h^*^ = {R^^^ / R^^^)"^ ■ bp^^ . The associated toroidal 
magnetic field is equal to {Q.^^^ ■ R/V^^^) ■ B^^^ ■ sin 9 in the RSG wind, 
and to {QP^^ ■ R/V^^^) ■ B^^° ■ sin 9 in the BSG wind. It may be notable 
here that the magnetic field is virtually toroidal, as can be seen from the 
above. The poloidal field is truly minor and follows the meridional flow, 
which is, however, affected by the toroidal field distribution. 

Under these conditions, the ram pressure pv"^ , the thermal pressure P, and 
the magnetic pressure i?^/2^o at = 0.1 ly are 5.6 x 10"^'^, 2.3 x 10"^'^, and 
7.1 X 10^-'^'^ respectively, expressed in units of N ■ for the RSG wind, 
and 9.3 x 10"^, 3.3 x 10-^°, and 2.3 x IQ-^i, respectively for the BSG wind. 

There are 960 x 960 grid points for physical domain 1 ly x 1 ly. For the 
initial condition, the isotropic RSG wind, the split-monopole poloidal and 
the associated toroidal magnetic fields are set in a whole simulation box, 
except on the inner boundary, which is located at Rib = 0.1 ly from the 
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SN. The isotropic BSG wind and the BSG poloidal and toroidal magnetic 
fields are prescribed on the inner boundary. On this boundary, the BSG 
poloidal field is the same as the RSG one, while the BSG toroidal field is 
prescribed by the relation, {Q^^^ ■ Rib/V^^'^) ■ B^^'^ ■ sin 6. Although 
this toroidal field is smaller than the RSG one, some strong toroidal field 
is expected to be accumulated at the contact surface between the RSG and 
BSG winds. The time evolution of the interaction between the two winds 
is computed step by step using the two-step Lax-WendrofF scheme. 

3. Computational Results and Conclusions 

Figures 1 and 2 show the global patterns of the density and magnetic field 
intensity, respectively, when the forward shock propagates up to « 0.65 ly. 
A finite-density enhancement along the equator due to the equatoward mag- 
netic pressure force is found in the RSG wind in this figure (and also in figure 
3) . This enhancement is a natural consequence due to a magnetic-pressure 
effect, which contrasts to some enhancements under an ad hoc assumption 
in HD schemes by Luo and McCray (1991) and Wang and Mazzali (1992). 
It is worth noting that the equatorial enhancement is not found in the BSG 
wind, where the toroidal magnetic field is relatively weak. A dramatic 
magnetic-pressure effect is found at around the contact surface where the 
BSG wind sweeps up the RSG wind. That the density increases at the 
shell-like structured region of the interaction between the contact surface 
and the forward shock over all latitudes has been confirmed. The special 
sharp enhancement, which corresponds to the observed ring-like structure, 
is found on the equator. Its peak density is 3.2 x 10^ cm~^. 

Here, the role of the magnetic pressure is emphasised; if this pressure is 
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absent the accumulated RSG wind matter is not confined in a ring-like stru- 
cuture. Figure 2 shows that due to the sweeping-up effect of the BSG wind 
the magnetic field is also accumulated at around the contact surface. It 
is remarkable that the magnetic intensity decreases drammatically to zero 
at the equatorial edge of the shell, where the sharp density-enhancement 
is found in figure 1. This means that the surrounding magnetic pressure 
pushes and confines the RSG material to the equator from both the north- 
ern and southern sides, and thus the ring-like structure is formed. Between 
the reverse shock located at ~ 0.27 ly and the contact surface near to the 
forward shock we find a hill-like widespread magnetic field accumulation of 
the shocked BSG toroidal field. As is shown below, the magnetic pressure 
here also works on the BSG wind. 

In this paper our attention has been focused on the magnetic effect around 
the equatorial plane. Here, we note that, in addition to the equatoward 
force in low latitudes, there is also a poleward force of the magnetic pres- 
sure in the middle and high latitudes which results in a density enhancement 
around the polar axis (Washimi, Shibata 1993). Because the latitudinal de- 
pendence of the toroidal field is B^f, oc sin 6 for 0° < 9 < 9q in the initially 
prescribed RSG wind, the magnetic pressure also pushes the RSG wind 
poleward in the middle and high latitudes. In fact, due to this poleward 
pressure force along the shell- like enhancement, the density increases with 
latitude in the middle latitudes, but it decreases again close to the polar 
region, and thus a local density maximum of ~ 1.6 x 10^ cm"^ is found at 
9 ~ 15°. This is probably due to a dynamical effect of the BSG wind (see 
figure 3), which prevents a density accumulation just on the pole. 
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4. Discussion 

In this MHD simulation analysis we have obtained an enhanced density 
of 3.2 X 10^ cm~^ at the ring. This peak value will probably be revised 
to a greater one by a more sophisitcated analysis, and will be consistent 
with the observed density (2-3) x 10^ cm~^ (Lundqvist, Pransson 1991). 
A suppression of the peak density value could probably be caused by the 
numerical diffusion effect, although we have attempted to reduce this effect 
by using a large mesh number. In addition, if our calculation could be 
performed from very near to the progenitor, where the poloidal field would 
be a dipole-like, rather than a split monopolc-like configuration, the equa- 
torward concentration of the material in the RSG wind would occur in a 
much greater colatitudinal width (e.g., 45°-135°), because the amplitude of 
the associated toroidal field (oc sin 9 ■ cos 9) is maximum at 9 = 45° and 
135°. For that case, the equatorial enhancement would be greater than our 
computational value, for which the split-monopole of the angle 85°-95° is 
conventionally used. In addition, if some cooling effect could be taken into 
account in the simulation, the density enhancement could be larger. 

We note that the idea of a magnetic field effect is consistent with the 
recent radio observation of a supernova remnant, SNR 1987A, detected by 
Staveley- Smith ct al. (1992) at about 1200 days after the explosion. This 
radio emission is explained by the collision of the supernova blast wave with 
the shocked blue wind (reverse shock at about 0.27 ly where some magnetic 
field increase and a small density change are found in figures 2 and 1, re- 
spectively). This position corresponds to the averaged expansion speed of 
the supernova ejecta ~ 0.08 ly which is consistent with the estimation by 
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Shigeyama and Nomoto (1990). The estimated magnetic-energy density by 
the minimum-energy argument is ~ 4 x 10~^/~^/^ A'^ m~^, where / is the 
fractional volume of the radiating acceleration region, suggesting a magnetic 
field of a few milli-Gauss or more (Chevalier 1992; Ball, Kirk 1992). This 
field intensity is consistent with an intensity of 2 • 10~^ Gauss obtained 
between the reverse shock and the contact surface shown in figure 2, if one 
takes into account a further enhancement of the field due to the sweeping-up 
process by the supernova blast wave. When the SN ejector collides with the 
ring at the end of this century or at the beginning of the next one, we can 
also expect more intense radio emission at rather middle and high latitudes 
where the magnetic intensity is greater, rather than at the equator where 
the ring-like structure is located. 

In addition, we note that the expanding speed of the ring in our present 
simulation is 32 km which is higher than the estimated speed of 10.3 km s 
(Crotts, Heathcote 1991). The resulting speed is a natural consequence of 
our parameters; we have used some conventional speeds of the RSG and 
BSG winds, 10 and 550 km respectively. We guess that these values 
should be lower to explain the ring speed. 

Chevalier and Luo (1994) studied the magnetic effect in the formation 
of the ring. Their idea is that the magnetic pressure and tension deform 
the shell to give an hourglass shape and the ring. However, their model 
is too simplified. Since both the latitude-dependence of the toroidal field 
strength, which has been shown to be essential for ring formation, and the 
inertial effect are not taken into account, their dynamical process for the 
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formation of the ring-like structure is much different from ours. 

It is interesting to note that the local-density maximum around the pole 
axis in our simulation constitutes additional rings in both the northern and 
southern hemispheres, which can corespond to the newly observed pair of 
rings. The opening angle ~ 15° and the diameter ~ 0.15 ly are different 
from the observed values (~ 45°, ~ 1.4 ly); however, these quantitative dis- 
crepancies could be overcome in future analyses if some conditions, such as 
high-speed wind flow in high latitudes, as observed in the solar wind, are 
additionally taken into account. 
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Figure Captions 

Fig. 1. Global density pattern for the region outside R = 0.1 ly (inner 
boundary) when the forward shock reaches 0.65 ly. The shell-like 
structure of the enhanced density is found, as well as the sharp peak 
on the equator corresponds to the ring-like structure. Its peak value is 
!^ 3,200 cm~^, while the peak value at 6' pa 15° is «i 1,600 cm~^. The 
positions of the reverse shock and the inner boundary are indicated by 
'RS' and 'IB', respectively. 

Fig. 2. Global toroidal magnetic field intensity pattern for the region outside 
R = 0.1 ly (inner boundary) when the forward shock reaches ~ 0.65 ly. 
Note the great field accumulation in the shocked RSG wind. The sharp 
decrease in the field accumulation appears at the equator, just where the 
sharp density peak in figure 1 is located. A hill-like accumulation of the 
magnetic field is also found between the contact surface and the reverse 
shock. The positions of the reverse shock and the inner boundary are 
indicated by 'RS' and 'IB', respectively. The plot does not show the 
exact equatorial plane where the toroidal field vanishes, but shows down 
to the grids just above the equator. 

Fig. 3. Equi-density pattern (solid curves) and wind velocity (arrows). 
The positions of the forward and reverse shocks and the inner boundary 
are indicated by 'FS', 'RS' and 'IB', respectively. The highest contour 
level is 10^-^cm~^, and the contour levels go down with the step 0.1 in 
logarithmic scale. 
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